Oersted fields and current density profiles in spin-torque driven magnetization 
dynamics — Finite element modelling of realistic geometries 
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The classical impact of electrical currents on magnetic nanostructures is analyzed with numerical 
calculations of current-density distributions and Oersted fields in typical contact geometries. For the 
Oersted field calculation, a hybrid finite element / boundary element method (FEM/BEM) technique 
is presented which can be applied to samples of arbitrary shape. Based on the FEM/BEM analysis, 
it is argued that reliable micromagnetic simulations on spin-tranfer-torque driven magnetization 
processes should include precise calculations of the Oersted field, particularly in the case of pillar 
contact geometries. Similarly, finite-element simulations demonstrate that numerical calculations of 
current-density distributions are required, e.g., in the case of magnetic strips with an indentation. 
Such strips are frequently used for the design of devices based on current-driven domain- wall motion. 
A dramatic increase of the current density is found at the apex of the notch, which is expected to 
strongly affect the magnetization processes in such strips. 

PACS numbers: 85.75.-d, 75.40.Mg, 41.20.-q, 41.20.Gz 



INTRODUCTION 



Current-driven magnetization processes have recently 
evolved to probably the most active field in magnetism. 
Recent discoveries demonstrating that, by exploiting the 
electron spin, the magnetization state of a nanomagnet 
could be influenced directly by electrical currents instead 
of external magnetic fields [J 0, HI, 0] have given rise to 
numerous extensive studies on the spin-transfer torque 
(STT) effect. Magnetization dynamics induced by STT 
@, Q is attractive for both technological aspects and 
fundamental physics. Unlike the field-induced magnetic 
switching, which in the case of nanodevices is connected 
with the technological difficulty of focussing magnetic 
fields on very small length scales, the STT effect pro- 
vides the possibility to switch, e.g., individual magnetic 
nanoelements in integrated circuits or in densely packed 
arrays Q. Moreover, the current-induced magnetization 
dynamics differs qualitatively from the field-induced dy- 
namics. A striking example thereof is the excitation of 
high-frequency magnetic oscillations induced by means 
of DC spin-polarized currents [l|, 0; an effect without 
analogy in the field-driven dynamics. The further pos- 
sibility of tuning the frequency of these oscillations by 
varying the current strength bears a high potential for 
future applications. 

The STT effect and the related fascinating features of 
current-induced magnetization dynamics have been ex- 
tensively studied by various groups over the last years. 
Micromagnetic simulations have nowadays developed to 
such a high level of accuracy and reliability that they 
provide an important tool for the design and for the un- 
derstanding of magnetic nanostructures and their mag- 
netic properties. While the STT effect has been implc- 
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mented in several micromagnetic codes @, [H, , there is 
often still a quantitative discrepancy between simulated 
and experimental results in the case of current-induced 
magnetization dynamics [j"lT | . This differs from the static 
[l2l [H[ and the field-driven cases [3| , where countless 
examples can be found showing excellent agreement be- 
tween experiment and simulation, thereby demonstrating 
the predictive power of classical micromagnetic simula- 
tions. 

One of the possible sources for such disagreement is 
connected with the Oersted field. Current-induced mag- 
netization processes require high current densities, typi- 
cally larger than about 10 11 A/m 2 . Such current densi- 
ties are inevitably connected with sizable magnetic fields, 
known as Oersted or Ampere fields [l5|. In addition to 
the usual external and internal effective fields entering 
the Landau-Lifshitz-Gilbert equation of motion, and in 
addition to the STT term, these Oersted field also af- 
fect the magnetization. The importance of the Oersted 
field in spin-transfer induced magnetic switching process 
has recently been pointed out by various groups [l6l [l7| . 
In view of the aforementioned high accuracy of classi- 
cal micromagnetic simulations, a numerical method is 
required in order to calculate the Oersted field contri- 
bution with high precision when STT effects are studied. 
In this article a method based on a combination of the 
finite-element method (FEM) and the boundary element 
method (BEM) is presented which allows for fast and pre- 
cise calculations of the Oersted field in arbitrary contact 
geometries. The hybrid FEM/BEM scheme has several 
advantageous features, including the simple considera- 
tion of inter-particle interactions resulting from Oersted 
fields generated, e.g., in complex nanocircuits [l8j or ar- 
rays of current-carrying nanopillars [19( . 

The paper is structured as follows. First, the basic 
equations and the numerical method for the Oersted field 
calculation are peresented. Subsequently, in section HVl a 
method to calculate current-density distributions in non- 
trivial geometries is outlined. Finally, a few examples 
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on the application of these methods to structures that 
are currently used for intensive studies on current-driven 
magnetization processes are presented in section Ivl 



II. BASIC EQUATIONS 

Consider a current-carrying sample of known size and 
shape. Let us further assume that the current density 
distribution j(r) is known for every point r in the volume 
V of sample. The calculation of j(r) will be discussed 
in in section IIVI The Oersted field H c connected with 
a current density distribution j(r) can be obtained by 
direct integration over the volume of the current-carrying 
particle: 



H(r) = Vx -L / T^rrdV' 
w 47r J \r - r'\ 

= Jj(r')x(V'G) dV 



where 



G(r - r') = 



1 



Air r 



is Green's function with 



AG 



-8(r 



(1) 
(2) 



(3) 



(4) 



and 6(r — r') is the Dirac delta function. Note that 
V'G = — VG. If the volume of the sample is discretized 
into finite elements, it is relatively simple to perform such 
a numerical integration. However, the direct integration 
according to eq. dJ) has a number of disadvantages from 
a numerical point of view. The volume integration over 
Green's function may give rise to serious problems of ac- 
curacy: Simple integration schemes can lead to inaccu- 
rate results because the value of Green's function can 
vary significantly within a single finite element contain- 
ing a point r' if the viewpoint r is close to that element. 
Moreover, the computational effort involved with this di- 
rect integration is very large: The calculation of each 
component of H c at a single point r requires an inte- 
gration over the whole volume V of the conductor. A 
more convenient approach in terms of accuracy and effi- 
ciency consists in numerically solving partial differential 
equations instead of the direct integration according to 
eq. flU). In the following, the basic steps for a FEM/BEM 
formulation of the Oersted field problem are outlined. 

The fundamental equation for the calculation of Oer- 
sted fields is Ampere's law: 



V x H = j 



(5) 



In a current-carrying ferromagnet with magnetization 
M, the magnetostatic field H is constituted by two 
parts: 



where H s is the demagnetizing field (also called dipolar 
field of stray field) created by the magnetic surface and 
volume charges, and the field H c is the magnetic Oersted 
field (also called Ampere field) due to the current den- 
sity j. The following equations apply to these two static 
magnetic field contributions: 



V x H s = 
VH, = -VM 



VxH c =j (7) 
VH C = (8) 



The calculation of H s is one of the central parts of 
any micromagnetic code. Several powerful methods have 
been discussed for the numerical calculation of the dipo- 
lar field [13, and this part of the problem can now 
be generally considered as solved. Due to its similarity 
with the method described later, one particular efficient 
method to calculate the demagnetizing field should be 
pointed out here, namely the hybrid FEM/BEM algo- 
rithm presented by Koehler and Fredkin [22|]. Hence- 
forth, only the contribution from the electric current H c 
shall be considerd, and the subscript "c" is omitted for 
simplicity. 

Equations and |(8]) yield 



AH = -Vxj 



(9) 



In the region outside the conductor, the source term is 
zero, i.e., 



AH = 



(10) 



Thus, each component of H satisfies an equation of the 
Poisson form inside the conductor and of the Laplace 
form on the outside. Equations ([9]) and (fT0|) describe an 
open boundary problem. Such problems are characterized 
by the absence of well-defined boundary conditions at the 
sample surface. The solution is uniquely defined by the 
condition of "regularity at infinity", i.e., a Dirichlet-type 
of boundary condition for r — > oo: 



lim H{r) = 



(11) 



H — H, + H, 



(6) 



A possible approach for the consideration of such 
boundary conditions consists in attempting to expand 
the computational region to "infinity", e.g. by applying 
bijective transformations to map the infinite volume sur- 
rounding the sample onto a volume of finite size, which 
can then be discretized with finite elements 23]. Such 
transformation methods can suffer from accuracy prob- 
lems because the discretized spatial transform effectively 
corresponds to a truncation of the computational region 
at more or less large distances. Moreover, these methods 
require the external region to be free of charges, thereby 
precluding the possibility of calculating the field in the 
case of interacting current carriers. An alternative ap- 
proach, which solves these problems, is the use of a hy- 
brid finite element / boundary element method. In the 
boundary element method, the fundamental solution of 
a partial differential equation is used. In this case, the 
fundamental solution is given by Green's function, which 
automatically fulfils the boundary condition at infinity. 
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III. HYBRID FEM/BEM FORMULATION 

The first step in developing the FEM/BEM formula- 
tion consists in an analysis of the properties of the solu- 
tion at the sample boundary. From eq. ^ a jump con- 
dition for the normal derivative of H c at the boundary 
of the current carrier is easy to derive: 



n( Vff x | in - VH x \ out ) = n(e x x j) 



(12) 



This equation correspondingly holds also for the y and 
z components of H . In eq. (|12p e x is the unit vector 
along the x axis and n is the surface unit vector directed 
outwards. Note that this condition represents a discon- 
tinuity of the gradient of H x at the boundary. However, 
it does not provide any information on the value of the 
inward limit of the gradient of H x . This condition can 
therefore not be used directly as a Neumann boundary 
condition at the boundary d£l as it would be required 
for a unique solution of the differential equation in the 
region Q. 

In order to obtain useful boundary conditions at the 
sample surface, the Oersted field can be split in two parts, 
H = JlW + H (2 l These fields shall have the following 
properties: (i): Outside the current carrier, the part H^' 
is zero, (ii): The part H^ 1 ' satisfies Poisson's equation 



-e x ■ (V x j) 



(13) 



with Neumann boundary conditions 



n ( Vff« 



= n (e x x j) (14) 

(iii): The part satisfies the Laplace equation 

Ai4 2) = (15) 

and its derivatives are continuous along the boundary. 

The boundary conditions required for the solution of 
eq. (fl~5)) will be discussed later. The jump condition for 
the individual components H^ and H^ at the bound- 
ary 



HP 



H& 



(16) 



results from the condition that the field H must be con- 
tinuos at the surface. The splitting of the field H x in 
two parts is schematically shown in Fig. [1] The advan- 
tage of this splitting is that it is possible to extract useful 
boundary conditions required for the solution of Poisson's 
equation: By setting Hx to zero outside the sample, 
the discontinuity condition (fP2|) is converted into a Neu- 
mann boundary condition (|14[) that can be used to solve 
eq. (fT3")) . With the finite element method, it is relatively 
easy to solve the Poisson equation (fl"3)) with given Neu- 
mann boundary conditions ([HI) . After that, the some- 
what more complicate part needs to be addressed, i.e., 
finding the values of H 1 - 2 ' at the particle surface. These 



internal region 



external region 




H=0 



FIG. 1: (Color online) The field H is split in two parts. The 
angle a represents the known value of the discontinuity of the 
derivative of H at the boundary (eq. 1 1 2 p . The part Hi of the 
solution is chosen such that its gradient at the boundary cor- 
responds to a. Moreover, Hi is equal to zero in the external 
region. Both, Hi and H2 are generally discontinuous at the 
boundary. The sum, i.e., the field H is however continuous. 
The gradient of the part H2 is constant along the boundary, 
as sketched by the angle /3. 



values will provide the Dirichlet boundary conditions re- 
quired for the solution of Laplace's equation (fl~5)) . For 
this task, the boundary element method is used. 

Multiplying eq. (U|) with Hi and eq. (fTB"]) with G yields, 
after subtraction and integration over the volume V of 
the region f2: 

J GAHfdV - J HfAGdV 

= J H*5{r - r') dV - J e x (V x j) G dV 

(17) 

which, by virtue of Green's theorem, transforms into 

+ in(e,xj)GdS'-e x [jx(V'G)dV 

The surface integrals § dS 1 extend over the boundary 
dfl of the region f2, i.e., the sample surface. Inserting 
the boundary condition for Hi 



dHl 
dn 



n (e x x j) 



(18) 



the equation simplifies to 

-j> Hf ^ dS' = J HfS(r-r') dV'-e x J jx ( V'G) dV" 

(19) 
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The last term can be identified as the right-hand side of 
eq. (JTJ), yielding 

-j H ^dS' = J H?5(r-r')dV'-[HZ(r)+HZ(r)] 

(20) 

The integral J HfS(r — r') dV is trivial as long as the 
point r is located inside or outside the volume V. The 
situation requires more attention when r is a point on the 
boundary. In this case, the inward limit is taken by intro- 
ducing an infinitesimal distance of the point (located in 
the internal region) to the surface. The resulting integral 
over the delta function is in this case 



HfS(r-r')dV' = ¥-Hf(r) (21) 



where ^ is the solid angle subtended at the boundary 
point r. 

Hence, an equation is obtained with which H' 2 (r) can 
be calculated at each boundary point: 

H*(r) = j> H?(r')^ dS' + ( ^ - l) Hf (r) (22) 

Accurate numerical methods to perform this integral by 
means of BEM are discussed in Ref. l 24]. In principle, 
eq. (j2"0|) would be sufficient to calculate H 2 at any point 
inside (and outside) the volume, so that the resulting 
total field H = Hi + H 2 could be calculated at any dis- 
cretization point. However, the calculation of H 2 at one 
point r requires in this case an integration over the whole 
surface. This approach would thus have similar disadvan- 
tages as the direct integration according to eq. ([1}. By 
using equation (|22|) instead, we obtain at relatively low 
cost the boundary values of H 2 . Having these Dirichlct 
boundary conditions, it is easy to solve eq. (|15[) in the 
volume V by means of the FEM. 

It is noteworthy that eq. (|2"2"|) has the same form as it 
has been used in the FEM/BEM scheme described in 
Ref. [12] to calculate the scalar magnetic potential of 
a ferromagnet. From the viewpoint of implementation 
into a program, this means that the matrix required for 
the numerical calculation the values of H 2 according to 
eq. (|22p is already available in the micromagnetic code if 
the FEM/BEM scheme presented in Ref. [22| is used for 
the calculation of the magnetic scalar potential. 



IV. CALCULATION OF CURRENT DENSITY 
DISTRIBUTIONS 

Unless the geometry of the current-carrying sample is 
trivial, the current density distribution j(r) needs to be 
determined numerically prior to the calculation of the 
Oersted field. The starting point for the calculation of 
the current density distribution is Ohm's law 



E is the local electric field and the conductivity a is 
assumed to be a scalar. The electric field is the gradient 
field of the electrostatic potential U, so that E = —'VU. 
Charge conservation yields 'VE = and thus Vj = 0, 
which ultimately leads to 



V(aVU) = 



(24) 



This equation has the form of the stationary diffusion 
equation and converts into the Laplace equation for U if 
a is homogeneous. 

Elliptic differential equations of the type are rou- 
tinely solved numerically with FEM. However, appropri- 
ate boundary conditions must be specified to obtain a 
unique solution. In this case, the boundary conditions 
are given by the fact that the current is not flowing per- 
pendicular to the sample surface j ■ n = except for 
the leads, where a known current density is entering and 
leaving the sample. In the case of two contact leads of 
the same size, the boundary condition is j ■ n = ±jo at 
the leads, leading to 



dU 
dn 



=Fjo/<7 at the leads 







at the rest of the boundary 



(25) 



j = aE 



(23) 



In a more general case, only the value of the total cur- 
rent / flowing through the sample is known. The value 
of the current density at the positive and negative leads 
then needs to be determined by dividing the total cur- 
rent flowing through each lead by the area of the contact 
region. Note that eq. ([24]) with given Neumann bound- 
ary conditions is as easy to solve numerically for homoge- 
neous conductivity a as it is in the case of inhomogeneous 
conductivity cr(r). 

The method described here for the calculation of cur- 
rent density distributions is generally well known and has 
been applied in various previous publications (25l . l26l [2^ | . 
It is included in this manuscript mainly for complete- 
ness, because such current density distribution calcula- 
tions usually are a necessary prerequisite for the calcula- 
tion of Oersted fields. 



V. EXAMPLES 

The methods outlined above can be applied to vari- 
ous problems which are of high importance for modern 
research topics in the field of nanomagnetism. In this 
section, a few examples are presented. These include the 
current density distribution in a thin strip with an in- 
dentation ("notch") and the calculation of the current 
density distribution and of the Oersted field in a pillar 
contact geometry as it is used to study high-frequency 
excitations in nanomagnets [l], [28| . 
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A. Current density distribution in a thin strip with 
a notch 

Ferromagnetic thin strips with width in the sub-/itm 
range and thicknesses of the order of a few ten nm have 
attracted much interest over the past years. One reason 
is the particular type of head-to-head domain walls that 
occur in such structures [2^ . Magnetic strips have been 
proposed to be used in future nanomagnetic devices, in 
which the head-to-head domain walls would serve as units 
of information that could be processed in logical devices 
[3fl l3l| . It has further been shown that the domain walls 
in such strips can be displaced by means of the STT 
effect if an electric current of sufficient strength is flowing 
through the sample [H, E| • Since the domain walls can be 
displaced continuously along a magnetic strip, it is of 
practical importance to gain control of their position, so 
that they are located and displaced between well-defined 
positions in the strips. This can be achieved by means 
of indentations or notches in the strips, which act as an 
attractive potential for the domain walls [HJ . Depending 
on their type, the domain walls are then either located 
exactly at the notch or in its close vicinity (32l . l33j . 

The influence of notches on the micromagnetic con- 
figuration has been studied intensively and reported in 
several publications. An equally important question for 
the study of current-driven domain wall motion in such 
indented strips is the influence of a notch on the current 
density distribution. 

To address this question, the current density distri- 
bution can be simulated with finite element modelling 
according to the method described in the previous sec- 
tion. The model used for the simulation is a 1.2 /im long 
strip (width: 100 nm, thickness: 20 nm) with an inden- 
tation in the center. The indentation reduces the width 
of the strip to 50 nm in the narrowest part. The angle of 
the indentation is 45°. The region of interest of the finite 
element mesh is shown in the inset of Fig. O In the vicin- 
ity of the notch, the tetrahedral mesh is locally refined 
to increase the numerical accuracy. A total current of 
1 = 2 mA is flowing along the wire, and a homogeneous 
electric conductivity is assumed. 

The topographical representation in Fig. [2] displays the 
local value of the current density. As expected, the cur- 
rent density increases at the constriction. Much more 
significant than the reduction of the width is, however, 
the impact of the notch. Where the wire is not con- 
strained, the current density is j = 10 12 A/m 2 . At 
the apex of the notch, the value increases drastically to 
j = 5.5 • 10 12 A/m 2 , while in the opposite, flat part of the 
strip it only increases up to j = 1.55 ■ 10 12 A/m 2 . Hence, 
in the narrow part of the strip, the current density distri- 
bution is highly inhomogencous. Over the small distance 
of 50 nm it changes almost by a factor of four. The val- 
ues of the local current density obviously scale linearly 
with the applied current. The value of 2 mA has only 
been chosen here as an example since the resulting cur- 
rent density values are of the order of those reported in 




FIG. 2: Current density distribution in a thin, flat strip with 
a notch. The inset shows a part of the finite-element mesh 
used for the analysis close to the constriction. The main im- 
age shows a perspective view on the region of interest of the 
indented strip and a topographical representation of the mag- 
nitude of the local current density. The strip is contacted such 
that at large distances from the notch the current flows par- 
allel to the strip. The local current density is dramatically in- 
creased near the apex of the notch. Contratry to this, the cur- 
rent density on the flat side of the constriction is only slightly 
elevated compared with the value in the unconstrained parts 
of the strip. 

corresponding experimental studies. The current density 
profile is independent of the value of the applied current. 
Hence, the inhomogeneities of the current density distri- 
bution are directly connected with the sample geometry, 
and not with the value of the applied current. The profile 
is moreover invariant with respect to scaling. 

The knowledge of this drastically inhomogeneous cur- 
rent density distribution at a notch is expected to be 
of utmost importance for the design and for the under- 
standing of devices based on current-driven domain-wall 
displacement. The huge differences of the local current 
density in the constriction region are expected to have a 
significant impact on the local strength of the STT ef- 
fect, the onset of electromigration processes and on the 
heating of the sample. The seemingly plausible approx- 
imation of homogeneous current density distribution in 
the constriction region hence clearly appears to be inad- 
equate, and should be dropped. 

B. Oersted field calculation - comparison with 
analytics 

The method presented in section IIIII to calculate the 
Oersted field for a given current density distribution can 
be tested by comparing the numerical results with known 
analytic solutions. A simple example is the magnetic field 
of an infinitely long current-carrying cylinder with homo- 
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FIG. 3: (Color online) Oersted field of a long, current-carrying 
cylinder with homogeneous current density. The current is 
flowing parallel to the wire axis. In this case, the field H only 
has an azimuthal component. The computed values (dots) 
perfectly match with the analytic solution (line). The grey 
area denotes the region inside the wire. 



geneous current density. The current is flowing parallel 
to the symmetry axis. If a wire with radius R is parallel 
to the z direction, with the current / flowing along e 2 
and (j) and p being the azimuthal and radial coordinate, 
respectively, the Oersted field is 

u r f P/R 2 (P < R) 

B(r) = poH(p) e(P = { (26) 

2n ( I/P (P > R) 

In Fig. [3] this analytical result (solid line) is com- 
pared with the computed values (dots) resulting from 
the FEM/BEM simulation in the case of an 8/Ltm long 
wire with R = 50 nm and homogeneous current densi- 
tiy j = 10 11 A/m 2 • e z . Excellent agreement is obtained. 
Numerical tests show that minor deviations can be fur- 
ther reduced by increasing the discretization density, but 
are not of practical importance. This example proves 
the correctness and accuracy of the hybrid FEM/BEM 
scheme. Obviously, the advantage of the FEM/BEM al- 
gorithm is its applicability to samples of complex shape, 
which cannot be calculated analytically. Such an exam- 
ple is given in the following subsection, where the cur- 
rent density distribution and the resulting Oersted field 
are calculated for a complex contact geometry as it is 
used to study current-induced stationary high-frequency 
excitations of nanomagnets. 




FIG. 4: Perspective view on the finite-element mesh used for 
the contact geometry. The ends of the 4/j.m long contacting 
strips, where the current is flowing into and out of the sample, 
are located outside of this frame, which is a magnified view 
on the region of interest where the pillar is located. The mesh 
is locally refined at the nanopillar. 



C. Current densities and Oersted fields in a 
nanopillar contact geometry 

A remarkable difference between the magnetization dy- 
namics induced by spin-polarized electrons as compared 
to the field induced dynamics is the possibility of gen- 
crating high-frequency, stationary oscillations in a nano- 
magnet with a DC current. This effect, which can occur 
when a sub-micron sized ferromagnetic thin-film element 
is exposed to a spin-polarized current flowing perpen- 
dicularly through its surface, has been predicted [f| Q 
theoretically and confirmed experimentally [l| . In exper- 
imental setups, the thin-film element is usually embed- 
ded into a pillar-shaped multilayer nanostructure, which 
is contacted by mesoscopic leads on the top and on the 
bottom, so that the electric current flows parallel to the 
pillar axis. Numerous experimental and numerical stud- 
ies on the current-induced magnetization dynamics in 
nanomagnets within a pillar contact have been reported 
recently (see, e.g., [34| and references therein). While 
the STT effect is driving the magnetization dynamics, it 
should be kept in mind that the Oersted field connected 
with the electric current may represent a non-negligible 
perturbation which could affect the magnetization dy- 
namics. Considering that also in this case the typical 
current densities are of the order of 10 12 A/m 2 , the Oer- 
sted field is expected to provide a sizable contribution 
even though the pillar dimensions are in the sub-micron 
range. 

To obtain precise values of the Oersted field in such 
complicated contact geometries, the approach consists in 
first simulating the current density distribution in the 
contact and subsequently calculating the magnetic field. 
Apart from the current flowing through the pillar, also 
the current in the leads (which act as what has become 
known as "strip lines" in the case of field-induced mag- 
netization dynamics) contributes to the total magnetic 
field acting on the nanomagnet. Fig. 2] shows the finite- 
element mesh used to simulate the current and field dis- 
tribution in a pillar contact geometry. In this exam- 
ple, the pillar diameter is 50 nm and the pillar height 
is 120 nm. It is contacted by strips of 500 nm width and 



FIG. 5: Simulated current density distribution in the contact 
geometry. The arrows represent the direction of the local 
current, and their length is proportional to the local value of 
the current density. For better visibility, only a selection of 
the computed data is shown. For the same reason, the current 
density distribution inside the pillar is not shown, since the 
much larger current density there would result in very long 
arrows in this representation. The top electrode is displayed 
as semi-transparent in order to visualize the current density 
distribution inside the leads. 



FIG. 6: Computed Oersted field resulting from the current 
distribution shown in Fig. [S] The circulation of the mag- 
netic field around the current-carrying regions is clearly visi- 
ble. The black arrows have a component directed outside of 
the sample, while the pale gray arrows are pointing towards 
the inside. Only the magnetic field at the surfaces of the leads 
is displayed. The magnetic field of the pillar is not shown for 
better visibility of the overall structure of the field. 



100 nm thickness. Each contacting strip is 4 /mi long and 
the top and bottom electrodes form an angle of 90° . The 
geometry roughly corresponds to typical experimental se- 
tups [23|, even though the contacting strips are often 
wider than shown here. The purpose of this simulation 
is rather to provide an example for the possibility of cal- 
culating current density distributions and Oersted fields 
in realistic geometries than to reproduce all details of a 
specific experimental setup. 

A current of 2 mA is flowing through the sample along 
the contacting strips. Since the total value of the current 
and the diameter of the pillar are known, the current 
density inside the pillar can be trivially calculated with- 
out any simulation (yielding j ~ 10 12 A/m 2 ). The same 
holds for the current density in the contacting strips at 
regions sufficiently far away from the contact (j = 4x 10 10 
A/m 2 ). The more interesting situation occurs in the 
leads in the vicinity of the contact to the pillar. The com- 
plex current density distribution in this region is shown 
in Fig. O 

The discretized form of the computed three- 
dimensional vector field j(r) can now be used as an in- 
put for the Oersted field calculation. By applying the 
FEM/BEM algorithm we then obtain the Oersted field 
in the whole contact geometry. The result is displayed in 
Fig. [6l where the field circulation around the contacting 
strips and, in particular, around the pillar can clearly be 
seen. For clarity, only the field in the contacting strips is 
displayed. 

For such experiments and their correct interpretation, 



the field distribution inside the pillar is decisive. The 
profile of the field inside the pillar is very similar to the 
one displayed in Fig. [3] The dominant component of the 
field is azimuthal, with a magnitude that inside the pillar 
increases linearly with the distance from the central axis. 
The peak value at the boundary of the pillar is in this 
case 28 mT, which certainly represents a non-negligible 
field. Compared with the field strength due to the cur- 
rent flowing in the pillar, the effect due to the current in 
the contacting strips can be considered as a perturbation, 
which provides only an asymmetry along a direction at 
45° with respect to the strips. Due to this asymmetry, 
the magnitude of the maximum and the minimum value 
of both, the x and the y component of the Oersted field 
differ by about 12%. Whether the influence of the cur- 
rent through the leads can be neglected depends on the 
specific experimental setup and the problem that is being 
investigated. While the field due to the contacting strips 
may in several cases be safely neglected, it is of essential 
importance to consider the precise height of the pillar, 
since it has a direct and decisive impact on the strength 
of the Oersted field. Reliable calculations of Oersted 
fields in nanopillar geometries can only be obtained if 
the height of the current-carrying pillar is known. More 
details on the influence of the pillar height on the Oer- 
sted field and its impact on current-driven magnetization 
dynamics will be discussed in a forthcoming article. 



8 



VI. CONCLUSION 

The classical interaction of electric currents with the 
magnetization of ferromagnets is given by the Oersted 
field. Its importance should not be underestimated 
in theoretical studies on the emerging field of current- 
induced dynamic magnetization processes driven by the 
STT effect. As shown in this study, the values of the 
Oersted field obtained by using finite-clement models of 
typical experimental contact geometries can be as high 
as a few tens of mT. In view of these large fields, it ap- 
pears likely that the details of several current-induced dy- 
namic magnetization processes that are currently studied 
by many groups are the result of a combination of both, 
the STT effect and the Oersted field. Thus, for a reliable 
interpretation of experimental data by means of micro- 
magnetic simulations, the Oersted field may a priori not 
be neglected. The hybrid FEM/BEM method presented 
in this article provides a flexible and accurate tool to 
calculate the Oersted field, and thus to clarify its impor- 
tance and its influence on the dynamics of nanomagnets 
in STT studies. The value of the Oersted field depends 
sensitively on the geometry of the contact setup. For in- 
stance, in pillar geometries, the pillar height is of decisive 
importance. 



Also in the case of current-induced domain wall dis- 
placement in ferromagnetic strips, some seemingly plau- 
sible simplifying assumptions may not be applicable. Mi- 
cromagnetic simulations should in this case be extended 
by a routine to calculate the current-density distribution 
like, e.g., the FEM scheme described in this article. On 
the example of the frequently used case of indented thin 
strips, the simulations have shown that geometric varia- 
tions can have a dramatic impact on the current-density 
distribution in thin strips. Such a notch can give rise to 
a strongly inhomogeneous current-density profile along 
the strip width connected with a strong, localized in- 
crease of the current density. It is therefore not sufficient 
to assume that the increase of current density in such an 
indented strip simply correlates with the reduction of the 
cross-section at the notch. 
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